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SUMMARY 


The SFD (sudden frequency deviation) events of May 18 and 19, 
1973, due to the solar flares as observed from the NASA/Marshall 
Space Flight Center high-frequency doppler sounder array system In 
Huntsville, Alabama, are Investigated. The results are compared with 
those observed at Table Mountain near Boulder, Colorado and at the 
University of Hawaii. 
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I. INTRODUCTION 


A sudden frequency deviation (SFD) is a type of sudden ionospheric 
disturbance (SID) causeo by bursts of X-rays and extreme ultraviolet (EUV) 
radiation from solar flares. It consists of a rapid change in the received 
frequency of a high frequency radiowave propagated from an ultra-stable 
transmitter and reflected at altitudes in the 120-300 Km range or higher 
in the ionosphere. 

The characteristic frequency derivation results from a time-varying 
ionospheric electron density produced by the time-varying ionizing radiation 
from a solar flare and is only sensitive to the impulsiveness of the 
flare radiation. Conversely, one can deduce the flare enhancement of the 
ionizing radiation from the sudden frequency deviation event. 

Thi6 report presents the results of the study of the SFD events of 
the flares of May 18 and 19, 1973, as observed from NASA/MSFC high- 
frequency doppler sounder array system in Huntsville, Alabama. The 
results are compared with those observed at Table Mountain near Boulder, 
Colorado and at the University of Hawaii. 



II. FLARE INDUCED EFFECTS IN THE IONOSPHERE 


II. 1 Solar Flare 

Solar flare generally means a burst of radiation occurring In the 
chromosphere of the sun. A flare Is customarily observed In the Ha light 

O 

(6563 A) as well as X-Ray emissions, although some flares are seen In white 
light also. The development of a flare Is usually as follows: a rapid rise 
to peak Intensity, a brief period of peak intensity followed by a steady 
decline. During a flare, there Is a temporary brightening of a small 
port i on of the solar chromosphere, usually near a sunspot. The active 
areas are associated with sunspots, plages, filaments and high magnetic 
field gradients. 

As expected, the flare is a complex phenomenon comprising diversified 
dynamical processes occurring in the solar atmosphere. Flares are frequent 
occurrences, particularly during the peak of the sunspot cycle. 

During most flares, there is a marked increase in the flux of solar 
ionizing radiation in the far ultraviolet and soft X-ray regions of the 
radiation spectrum. There are also radio noise emissions and explosive 
outbursts of corpuscular emission ranging from relativistic to lower 
speeds. These radiations mark their imprints in the appearance of various 
ionospheric effects, such as enhancement in the electron densities in the 
D and E regions, auroras, magnetic storms, etc. 

Nearly all important extensive disturbances in the ionosphere 
are associated in some way or other with a flare on the sun. These 
disturbances are important from the point of view of radio communications 
because they often result in interruptions of communications. The 
absorption in the D-region is enhanced so much that intelligible radio 
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communication becomes impossible for periods lasting from a few minutes 
to several days. Also, the critical frequencies of the F2 layers are 
sometimes depressed during an ionospheric storm resulting in loss of 
the signal due to MUF (Maximum Usable Frequency) failure. The study 
of these disturbances is important for a better understanding of the 
interaction of the solar radiation with the earth's atmosphere and 
ultimately a better understanding of the solar flare Itself. 

The flare-induced effects can be broadly classified into two 
categories: (1) the simultaneous effects and (2) the delayed effects. 

The former effects are due to electromagnetic radiation enhancements 
and the latter due to solar cosmic ray particles, both relativistic and 
non-relativistic. The relativistic particles (mostly protons) having 
transit times from 15 minutes to several hours Induce polar cap absorption 
and cosmic ray enhancements. Slower particles (both protons and electrons) 
which take 20 to 40 hours to reach the earth are responsible for ionospheric 
storms, magnetic storms and visible displays of aurora. 

The fact that many of the ionospheric effects begin at the time 
of visual sighting of the flare indicates that they are due to electro- 
magnetic rather than corpuscular radiation. The occurrence of a flare is 
accompanied by the emission of radio waves, UV - and X-rays, all of which 
arrive at the same time with visible light. The various types of radio 
noise are thought to be produced by bremsstrahlung synchrotron emission 
and plasma oscillations in the solar atmosphere (Wild, 1964? Smerd, 1964). 

As a result of a solar flare. all layers of the ionosphere are 
affected, although the E layer is not strongly affected except for the 
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occurrence of sparadic E in certain ionospheric storms. The D and F 
layers are subject to much stronger effects. 

The general classification "sudden ionospheric disturbance" or 
SID includes the phenomena which accompany solar flares occurring 
simultaneously with the optical Ha flare (Dellinger, 1973). SID 
observations form two groups. The first group consists of SID effects 
observed in the ionosphere at D region heights below 100 Km caused by 
soft X-rays (10A) and Includes SWF (Short Wave Fadeout) , SPA (Sudden 
Phase Anomaly) , SCNA (Sudden Cosmic Noise Absorption) , SEA (Sudden 
Enhancement of Atmospherics) and so on. These SIDs provide useful infor- 
mation for studying the transient response of the poorly understood D 
region (Mitra, 1972). 

II. 2 D Region Effects 

The D region effects tend to involve large percentage Increases 
in ionization relative to preflare values. Satellite measurements have 
shown that the total solar X-ray emission can easily rise by several 

O 

orders of magnitude in the wavelength 2-10 A in a class 2 flare (Culhane, 
et. al., 1964). An increase in electron concentration in the D region 
by a factor of 10 was measured by the certlcal reflection technique 
(Belrose and Centiner, 1962). This clearly identifies the X-ray 
enhancement as responsible for the D region SID effects. The threshold 
effect results from cosmic rays and the solar hydrogen Lyman a line, 
the dominant source of ionization in the D region, until the soft X-ray 
flux raises this threshold value. 

The production of ionization in the D region gives rise to associated 
phenomena such as the short wave fadeout (SWF), sudden cosmic noise 
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absorption (SCNA) received on frequencies above the F2 layer critical 
frequency, the sudden phase anomaly (SPA) and the sudden enhancement 
of atmospherics (SEA) on very low frequencies. In the following, we 
shall give a brief description for these terminologies. 

a. Short Wave Fadeout (SWF) - At times communications on high 
frequency by skywave propagation over the daylight hemisphere of the 
earth are blacked out or strongly absorbed by enhanced ionization in 

the D region. There is a sudden decrease in the received signal strength 
beyond the normal fading of the HF radio wave reflected from the ionosphere. 
The association between these short wave fadeouts and solar flares was 
discovered by Dellinger (1935) . 

b. Sudden Cosmic Noise Absorption (SCNA) - During a flare, there is 
a sudden decrease in galactic cosmic noise received at frequencies high 
enough to penetrate the ionosphere yet low enough to encounter measurable 
absorption (Jansky, 1937). A fairly large flare may cause a 2 db increase 
in absorption of frequencies of 18 to 25 MHz which pass completely through 
the ionosphere without reflection (Shaln and Mitre, 1954). 

c. S udden Phase Anomaly (SPA) - The Increase in ionization produces 
a phase advance of low frequency radio waves reflected from the iono- 
sphere due to a lowering of several kilometers in the height of reflection 
(Budden and Ratclifte, 1937). 

d. Sudden Enhancement of Atmospherics (SEA) - During a flare, there 
is a sudden increase in signal strength of noise (atmospherics) from 
lightening storms around the world. These are recorded between 20 to 40k Hz 
(Bureau, 1937). Like the SPA, this indicates a change in very low 
frequency propagation conditions. 
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e. Solar Flare Effect (SFE) or Magnetic Crochet - SWFs are sometimes 
accompanied by transient variations in the earth’s magnetic field. The 
horizontal component of the geomagnetic field undergoes a sudden small 
change (McNish, 1937). This indicates the existence of electric currents 
in the D region ionosphere, due to increased ionospheric conductivity. 

II. 3 E and F Region Effects 

The solar hydrogen Lyman a line is the dominant source of the 
SID effects at heights above 100 Km in the well-understood day-time E 

O 

and FI regions caused by solar radiation in the 10 to 1030 A wave- 
length range. These SID effects involve much smaller percentage increases 
in the electron density but much larger absolute increases in comparison to 
those in the D region. 

Most of the results on the EUV flare radiation are deduced from 
ground-based measurements based on one observing technique where the 
flare effect is called a sudden frequency deviation (For review, see 
Donnelly, 1970). In this experiment, the frequency of a very stable 
radio signal is observed after reflection from the E or FI layers. An 
increase in the electron content of the ionosphere during a solar flare 
brings about a decrease in the phase path of the signals and induces a 
small transient change In frequency of the observed signal. Although 
the extra D region ionization produces large absorption effects, it 
contributes very little to the phase path changes measured in an SFD. 

These changes are mostly due to extra ionization deposited in the E 
and FI regions by longer wavelength X-rays and EUV radiation (Kanellakos, 
et al., 1962). 
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III. SUDDEN FREQUENCY DEVIATIONS 


111.1 History of Sudden Frequency Deviations 

The sudden frequency deviations due to solar flares were Identified 
at a much later date than most other SIDs (Chan, et al., 1961). Chan and 
Vlllard (1963) used the term "sudden frr uency deviations" to describe 
the HF Doppler observations of SIDs. Davies, et al. (1962) derived the 
frequency deviation for the case of a thin enhancement region below 
the height of reflection and for the case where the enhancement of 
Ionization Is at the height of reflection. They showed that the frequency 
deviation Is Inversely and directly proportional to the transmitter 
frequency In the first and second cases, respectively. 

Statistical studies of SFDs have been made by Chan and Vlllard 
(1963) and Agy, Baker and Jones (1965). Recent investigations have been 
made by Donnelly (1970) and Donnelly, et al. (1974). The main assets of 
SFD observations are 1 sec time resolution, low cost, ground-based 
equipment and essentially continuous daylight coverage. 

111. 2 Characteristics and Wavelength Dependence of the SFDs 

In an SFD event, the received frequency of a high frequency radio 
wave, usually reflected from the F region of the Ionosphere, Increases 
suddenly, peaks and then decays to the transmitted frequency. The 
received frequency usually decreases below the preflare frequency from 
the peak and then increases back to the preflare frequency. At other 
times, no negative frequency deviation takes place. Some SFDs have 
more than one peak, and most SFDs have some fine structure. The 
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start-', ^-maximum time is typically about one minute and the peak 
deviation is usually less than 0.5 Hz. The size of the frequency 
deviation Is very small relative to the transmitted frequency, rarely 
more than one part per million. Unlike most other detectors of solar 
radiation, the SFDs are sensitive to the impulsiveness of the solar 
radiation and insensitive to the non-flare radiation. 

The SFDs differ from the other SID effects in that the frequency 
deviation is proportional to the time rate of change of electron density 
primarily in the E and FI regions produced by flare radiation in the 

O O 

10-1030 A wavelength range. Radiation enhancements in the 1-10 A range, 
although responsible for SFDs from the bottom of E layer and the sporadic 
E. contribute mainly to the D region electron density enhancements and 

O 

the resulting SID effects. Radiations at wavelengths greater than 1030 A 
are ineffective in ionizing any of the major constituents of the 
atmosphere, and are therefore, unable to contribute to any SID effect. 

O 

The sensitivity of SFDs to the 1-1030 A bursts depends upon the spectrum 

O 

of the l-j.030 A radiation, the solar zenith angle, time of the day and 
season, the preflare electron density as a function of height, the upper 
atmospheric neutral densities, etc. 

III. 3 Sensitivity of SFDs to the Solar Radiation 

The sensitivity S(X) of SFDs to incident radiation as a function 
of the wavelength X of the ionizing radiation is defined by 


S(A) 


Af (fy) 

R t Rjf A<f> (X) dX 
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where Af is the frequency deviation in Hz, A4> is the radiation enhancement 


in 


ergs cm -2 sec 1 A -1 , R^ is a dimensionless factor that accounts for 


the time-dependence of A<J>, R , accounts for the solar zenith angle 

A 

dependence and f^ is the equivalent vertical-incidence frequency of the 
SFD probing radio wave. 

Figure 1 gives a plot of S(A) for a radio wave reflected at 200 
Km altitude. The value of S(A) drops off rapidly with decreasing wave- 


length below about 1 A, because radiation at those wavelengths produces 
ionization mainly low in the D region, where most of the freed electrons 
are rapidly lost via attachment to form negative ions. 

Because of small ionospheric fluctuations that are always present, 

O 

the smallest 1-1030 A flux enhancement detectable by mid-latitude SFD 

observations under the best conditions is about 4 x 10 erg cm sec 

2 2 1 

and more typically 10 erg era - sec - . The frequency deviation is 

O 

linearly related to the 1-1030 A flux enhancement for enhancements less 
than 1 erg cm sec , above which a small non-linearity may occur after 
the EUV flux has been in progress for about 100 sec. 

Figure 2 illustrates the sensitivity of SFDs to the impulsiveness 

O 

of 1-1030 A radiation. The solid curve represents the tine dependence 


O 


of a hypothetical EUV burst of ionizing radiation in the 10-1030 A wave- 
length range, or the resultant rate of electron production with a linear 

rise to a peak in time l followed by an exponential decay with a 1/e 

m 

time constant t. set arbitrarily equal to tt t . The dashed curve 
d 2 m 

illustrates the time-dependence of the rate of change of electron density 
and hence the SFD for the case where the rise time t and the decay time 
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Figure 1 . SFD sensitivity to ionizing radiations as a function of radiation wavelength. 
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t , of the burst are much less than the effective electron-loss time 
d 

constant t In the E and FI regions of the ionosphere. In this case, 
c 

the SFD time dependence is essentially the same as that of the incident 
radiation enhancement, except late in the decay stage. The dash-dot 
curve illustrates the case when the rise and decay times of the burst 
of radiation are comparable to the effective electron-loss time constant. 
The SFD in this case has a time dependence similar to the EUV burst but 
with a marked negative deviation during the decay stage. Most SFDs fit 
into this case. The dotted curve illustrates the case when the burst 
time constants are much larger than the electron-loss time constant. 

Here the resulting SFD is smaller and flatter. The electron-loss time 
constant in the E and FI region is typically between 15 and 60 secs. 
Hence, SFD observations are relatively insensitive to radiation bursts 
with smooth rises of a few minutes or more. This illustrates the SFD 
dependence on the impulsiveness of the flare radiation. 
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Figure 2. The dependence of SFD's on the impulsiveness of the 
flare enhancements of ionizing radiation, i c * the 

effective electron-loss rate in the ionosphere. 
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IV. ANALYSES OF SFDs OF MAY 18 AND 19, 1973 

IV. 1 Description of Observed Events 

Figure 3 shows the map of the National Aeronautics and Space 
Administration/Marshall Space Flight Center high frequency doppler sounder 
system near Huntsville, Alabama, which was designed to study moving 
ionospheric irregularities including traveling waves excited by static 
ground tests of large rocket engines. Huntsville is located at 34.7°N 
latitude and 86.6°W longitude. 

The two SFD events as observed in Huntsville and analyzed in this 

report are shown in Figures 4 and 5. The event of May 18, 1973 started 

at 1525 UT (Universal Time) which corresponds to 9.25 a.m. local time 

and a solar zenith angle of 33.3°. The event showed two major peaks at 

1528 UT and 1530.7 UT and also negative frequency deviation between the 

two peaks. The maximum frequency deviation Af was .664 Hz at 1528 UT. 

max 

The prominent SFD event of May 19, 1973 started abruptly at 2243 UT 
corresponding to 4.43 P.m. local time in Huntsville and a solar zenith 
angle of 67.3°. The event had a single major peak at 2245 UT corresponding 
to Af ^ = 2.254 Hz. This large SFD event accompanied a IB Ha flare in 
McMath Plage Region 12352 at N09 E20 starting at 2233 UT and peaking at 
2245 UT (see Donnelly, et al., 1974). This event as well as the event 
of May 18 as observed in Huntsville appeared too smooth and devoid of fine 
structures. 

The same event (of May 19, 1973) observed near Boulder (40.2°N, 

105. 3 # W) at 137 MHz frequency and in Hawaii (21.3°N, 157. 8°W) at 10 MHz 
are shown in Figures 6 and 7 (Donnelly, 1976). This corresponds to a 
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Figure 3. Locations of transmitters, receivers, and ground paths for HF Doppler 

observations at NASA/Marshall Space Flight Center, near Huntsville, Alabama. 
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8. IS 1978 HUNTS VILLI ALABAMA MFSC PRBOUKNCY • 4.11 MHZ 
4B HZ PSAK PHIOUIHCY DWIATION At IBS7.BB UT 



UNIVERSAL TIME 


figure 4. Observed SFD event as a function of UT (Uni- 
versal Time) on May 18, 1973 at NASA/Marshall 
Space Flight Center high frequency doppler 
sounder system near Huntsville, Alabama (34.7* 
latitude and 86.6* W longitude). 
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NORMALIZED FREQUENCY DEVIATION 




9. 19 IOTS HUNTSVILLE ALABAMA MSFC FREQUENCY > All MHZ 
AT HZ PKAK FREQUENCY OBVIATION AT BS4B.II UT 



Figure 5. Observed SFD event as a function of UT (Uni- 
versal Time) on May 19, 1973 at NASA/Marshall 
Space Flight Center high frequency doppler 
sounder system near Huntsville, Alabama (34.7* 
latitude and 86. 6 # W longitude). 
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local time of 4:45 p.m. at Boulder and 12:45 p.m. In Hawaii, and there- 
fore, to smaller solar zenith angles. The maximum frequency deviations 
observed were also larger. The 4.0 MHz transmitter trace at Boulder was 
lost due to high flare-induced absorption and no comparison with that at 
Huntsville was possible. However, both the Boulder and the Hawaii data 
showed distinct fine structures. 

O 

IV. 2 Deduction of the Flare-Enhancement of the 10-1030 A Radiation Flux 

O 

The 10-1030 A flux enhancement A$ Is related to the flare-enhanced 
production of Ionization Aq(t) by 

A<f> (10-1030 A, t) * a sec x Aq(t) 

where a - 6.56 x 10“ 13 watt, m-sec. and x is the solar zenith angle at 
the midpoint of the propagation path at the time of the peak of the 
event (Donnelly, 1970). 

Aq(t) is given by 

dN AN 2 

» (t > - IT + ^ + B eff 

where x rr is the effective electron loss time-constant and the non- 
eff AN 

e 

linear term in B ec is small relative to . 

e£f T ef f 

O 

The procedure used to calculate the 10-1030 A flux enhancement 

for the results of this report was method three of Donnelly (1970). 

This method employs two main related simplifying assumptions: 

dN 

0 

(1) The time rate of change of electron density — — in the 

Ionosphere during the solar flare is approximately constant with height 

h » from just above the bottom of the E layer to the height of reflection h . 
o r 
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Figure 6. Observed SFD event as a function of UT (Uni- 
versal Time), on May 19, 1973 at Boulder, 
Colorado (40.2*N latitude and 105.3* W 
longitude) . 
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NORMALIZED FREQUENCY DEVIATION 


9.19 1971 MAUI TO U •• HAWAII PftIQUKNCY « 10.00 MHZ 
10.9 HZ PEAK PHIQUKNCY DCVIATION AT 1844.49 UT 



Figure 7. Observed SFD event as a function of UT (Uni- 
versal Time), on May 19, 1973 in Hawaii (21.3* 
N latitude and 157.8* W longitude). 
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(2) The electron loss time constant is approximately constant 

with height over the range from h^ to h^. 

The first assumption Is approximately satisfied at heights 110 to 

200 Km If the second assumption Is true and because the EUV flare spectrum 

« 

Is spread over the 10-1030 A range In continuum emission and In a large 

number of emission lines. The first assumption simplifies the relation 
dN e 

between Af and — j — (Agy, et al., 1965) 
at 

dN f c 

7? = if k. (h' - h o ) 

in the altitude range from h Q to h_, where k a 80.6 Hz 2 m 3 , f y is the 

equivalent vertical Incidence frequency, c Is the speed of light, and 

h f " is the vertical height of reflection. Transmission curves (Smith, 

1939) and lonograms taken just before the flare at a location near the 

HF Doppler propagation path are used to determine h^, h Q and f^ for 

each event and channel analysed. 

The values of T ,, and B ,, are determined from models of 
eff eff 

a (1) 8 (L) . The electron density as a function of height is 

eff eff 

calculated for each flare and HF Doppler observatory using the 
corresponding preflare lonograms and the methods of Wr.ght (1967). 

The true height of reflection h is then determined from N (h), where 

f 2 r e 
1 v 

N (h ) ■ -r— . Using several models of a _ _(h) or B(h) * 0 including 

e r * eff 

those >f Mitra and Banerjee (1971), T „„ and B are computed as a 

eff eff 

function of altitude from T(h) ■ [2a ,,(h) N (h) ] *, or l(h) ■ B _l 

eff e 

at the higher altitudes, and B(h) ■ a (h), or B(h) ■ 0 at altitudes 

611 

where 8 applies. Averages of i(h) and B(h) are computed over the height 
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range from h to h to determine T ,, and B ... 

or eff eff 

IV . 3 Results and Discussions 

O 

The normalized 10-1030 A flux enhancements as a function of time 
as deduced from the two SFD events by the method 3 of Donnelly (1970) 
are shown In Figures 8 and 9. The maximum Increase In the flux density 
max foi the May 18 event was 3.24 x 10“ 4 W/m 2 at 1533.7 UT. A$ 

max 

for the May 19 event was 8.79 x 10" 4 W/m 2 at 2245 UT. This compares with 

t 

1.56 x 10“ 3 W/m 2 as obtained rron. the Hawaii data (Donnelly, 1976). 

O 

The (10-1030 A) values computed using this method are estimated 
to be accurate to within a factor of four (Donnelly, 1970). The 

O 

Impulsive structure In A4> (10-1030 A) with time constants less than 

T ee are Insensitive to the height variations of T or to the a(h) and 
et t 

6(h) models, while the flux variations with time variations larger than 
are quite sensitive. Ionospheric variations unrelated to solar flare 
effects produce noise in the SFD observations that mainly affect the 
A$(10-1030 A) estimates during the slow flare emissions. S F D data 
provides good information on the impulsive emissions of flares, poor 
Information on the slow flare emission and no information on the non- 
flare solar radiation. 
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Figure 8 


UNIVERSAL TIME 

Best Estimate of 10-1030 A Flux Enhancements of 
1528 UT, May 18, 1973. 








DELTA T* (10- 1030 A) /DELTA PHI MAX (10-1030) 

RMALL T ~ I0-I030A SOLAR ENERGY FUiX ENHANCEMENT 



Figure 9. Best Estimate of 10-1030 A Flux Enhancements of 
2245 UT, May 19, 1973. 
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CHAPTER V. CONCLUDING REMARKS 


During this study, we investigated SFD events recorded on May 18 

and 19, 1973, from NASA/Marshall Space Flight Center high-frequency 

doppler sounder array system near Huntsville, Alabama. The May 19, 1973 

SFD event has been identified with a IB Ha flare in McMath plage region 

12352 at N09E20 of the solar surface. The maximum increase in the flux 

density has been deduced from these two observed SFD events, by using 

method 3 suggested by Donnelly (1970) . We found that the maximum 

4 * 2 

increase in the flux density for the May 18 event was 3.24 x 10 w/m 
at 2245 UT. This was a factor of - 2 to 5 smaller than the results 
obtained from the Hawaii data (Donnelly, 1976). Also, we have noted that 
the observed SFD events recorded in Huntsville are much smoother than 
those recorded at Boulder, Colorado and in Hawaii, which means that 
the observed results recorded in Huntsville, Alabama have less 
structure than those recorded in Boulder, Colorado and in Hawaii. This 
may be caused by local environmental conditions. However, the over- 
all features from the observation show no significant differences. 
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